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We present a method to study the time variation of the orbital parameters of a Post-Keplerian 
binary system undergoing a generic external perturbation. The method is the relativistic extension 
of the planetary Lagrangian equations. The theory only assumes the smallness of the external 
perturbation while relativistic effects are already included in the unperturbed problem. This is the 
major advantage of this novel approach over classical Lagrangian methods. 
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I. INTRODUCTION 



Since the early stages of classical celestial mechanics, a very large amount of efforts has been made in order to find 
exact or approximate solutions to the problem of N point -like interacting bodies. It is well known that a general 
solution exists only for N — 2 (Keplerian solution). In general, when N exceeds two searching for even approximate 
solutions becomes a very difficult task. Fortunately many N-body systems of interest in celestial mechanics may be 
considered as a 2-body problem, and the interaction of the other TV — 2 bodies can be regarded as a perturbation. 
In these cases the distance of the N — 2 bodies, together with the relative magnitude of masses, allows one to resort 
to standard perturbative methods (see for instance ch.6]). Moreover when a 2-body system suffers perturbations 
such as drag or radiation damping forces, oblateness of one of the two bodies and so on, the usual Keplerian solution 
must be viewed only as a zeroth-order approximation. 

In the framework of newtonian gravity, all these kinds of perturbations may be handled using the Lagrangian 
planetary equations through which the time dependence of orbital elements (otherwise constant) is achieved. This 
way the motion of the system is formally Keplerian but the orbital elements are allowed to vary with time. This 
procedure has been widely used from physicists and astronomers to study Newtonian binary systems; but what 
happens if relativistic effects are to be taken into account together with the external perturbation? A first way to 
approach the problem is to assume that relativistic effects and the external perturbation are roughly of the same order 
of magnitude; in this semiclassical standpoint they are both considered as a perturbation to Keplerian motion. 

We propose a novel procedure to account for relativistic effects in the unperturbed problem, in which the only 
perturbation is the external one. For this scheme we divide our work in two steps; first the relativistic two-body 
problem needs to be solved (this has already been done up to the IPN and 2PN orders beyond the classical limit by 
various authors (see, for instance, ||] and |^); this is the so-called post-keplerian solution); then our task, that is the 
main result of this paper, is to find out the relativistic version of Lagrangian planetary equations, giving the time 
dependence of relativistic orbital elements. 



II. POST-KEPLERIAN SOLUTION 



Let us consider a system of two bodies, moving under mutual gravitational attraction. Since we want to focus our 
attention on systems whose relativistic effects are not negligible, we have to introduce a Post-Newtonian parameter 
quantifying the relevance of these effects; this parameter is defined as follows 

def V 

epN = - (1) 

c 

In the above expression v is the typical speed of the bodies of the system and c is the speed of light; using relativistic 
Einstein equations it can be shown that the equations of motion up to the IPN order, i.e. retaining only the terms in 
(epAf)^, may be derived from the following Lagrangian, in a reference frame with the IPN centre of mass at rest (see 
for instance |^): 
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In the above formulae, ri and r2 are the radius vectors of the two bodies, mi and m2 denote their masses and is 
the position of the Post-Newtonian centre of mass (see jj, eqs.(5.4.5)-(5.4.9)]). In the following it will become clear 
the usefulness of the Hamiltonian approach with respect to the Lagrangian one; therefore we calculate the canonical 
momenta deriving from (|^): 

dCo V I v^v Gm - , Giym {v ■ x) 

The Hamiltonian function of the system is defined as usual: 

Ho=p-v-Co (4) 

Now, from (||) and we are in a position to express explicitly such an Hamiltonian in terms of the canonical variables 
X and p, obtaining: 
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Hence, Hamilton equations are: 
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The exact solution (called Post-Keplerian), up to IPN order, has been found by Damour and Deruelle (see and 
also for its extension to 2PN level). The Post-Keplerian solution may be written in a form very similar with 
respect to the classical one as follows: 
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In these equations cu, ft and i are Euler angles, defining the rotation that connects the observation reference frame 
with the intrinsic frame of the motion. In celestial mechanics they are usually referred to as argument of periastron 
(the angle in orbital plane from the line of nodes (see §4.4]) to the perihelion point), longitude of the ascending 
node (the angle measured from the positive x axis of the observer to the line of nodes) and inclination of the orbit 
(the angle between the orbital plane and the x-y plane of the observer), respectively. The other elements of the orbit 
are the semimajor axis of the ellipse a, the eccentricity e and the time of periastron passage T. 

III. IPN LAGRANGIAN BRACKETS 

In last section we have reviewed the post-newtonian solution of a binary system of bodies. Such a solution is 
quite correct provided that either the system is completely isolated or the external perturbation induces effects whose 
order of magnitude is less than or equal to the 2PN ones; if we quantify the weakness of external perturbation to the 
parameter eext, our last assumption is: 

eext < {epN^ (17) 

In the framework of the validity of the above expression the simple post-Keplerian solution is correct. When the 
external perturbation is so strong that the above expression does not hold true anymore, we have to take into account 
its effect on the IPN relative motion of the two bodies. To this purpose we have developed a perturbation method 
enabling us to calculate the time dependence of orbital parameters due to external perturbation. Such a method may 
be viewed as the relativistic extension of the planetary Lagrangian equations, which are widely being used in classical 
celestial mechanics to evaluate the time evolution of the orbital elements for Keplerian motion. Indeed it is our aim 
to calculate the time evolution of orbital parameters for Post-Keplerian motion. To summarize we assume that the 
solution of the perturbed problem has the same functional dependence on orbital parameters and on time explicitly 
as the unperturbed problem. 

In the presence of an external perturbation, the Hamiltonian describing the system can be written as: 

n = noix,p) + ni{x,p) (is) 

where Jio is the usual 2-body Hamiltonian, which is exact up to IPN order, while Hi describes the most general 
perturbation. This latter is a sort of disturbing function with the only difference that the usual classical disturbing 
function (see [Q, chs.5, 6] and Q if the external perturbation is due to the gravitational attraction of other bodies) is 
a perturbation to the Lagrangian, while the one concerned through this paper is a perturbation to the Hamiltonian. 
With the above assumptions, the solution fulfills the following equations: 

dx _ &Ho dp _ dHo 

dt ^ dp ' dt ~ dx ^ ' 

where x = x(a, e, T, w, 17, i, t) and p ~ p{a, e, T, uj, fl, i, t). 

As in the classical case, we need to find out the equations governing the time dependence of orbital parameters 
(a, e, T, w, f2, i) due to an external perturbation. First of all we have to calculate the relativistic Lagrangian brackets 
defined as follows & §9-4]: 
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where Cj with 1 < j < 6 are the orbital elements a, e, T, fi, i. On the strength of equations ( p9| ) it is easy to show 
that the Lagrangian brackets do not depend explicitly upon time (for the classical version of this property see 0, 
§6-2]). ^ ^ . . 

From definition ( 20|) and Hamilton equations of motions pertaining eq. (|l8| ) , we find, after straightforward calculation 
(see for instance p, §11-2]): 
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We remark that this latter equations are rather similar to the classical ones (see p. 133, eq.(6-19)]), where Hi plays 
the role of disturbing function (the sign is reversed since our IPN disturbing function is ascribed to an Hamiltonian). 

Since IPN Lagrangian brackets do not depend explicitly upon time, it suffices to evaluate them at the periastron 
(i.e. at t = T). In doing so we follow the same procedure as in the classical case (see for instance ]|l|, §6-3]), i.e. we 
expand in powers of t — T, up to second-order terms, the relative radius vector x and the canonical momentum p, 
then we compute Lagrangian brackets (po[); finally we set t = T. After very long but straightforward calculation we 
achieve the only non-vanishing brackets: 
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IV. IPN PLANETARY LAGRANGIAN EQUATIONS 

Using the above expressions for the brackets and inverting the system of equations ( ^l[) , we find the post-Newtonian 
Lagrangian planetary equations: 
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Once the external perturbation has been assigned, it is possible to evaluate its effects on the orbital parameters by 
solving the system (|2^)-(^); as in classical mechanics, the first-order time variation of orbital elements, is simply 
obtained after substitution of the unperturbed values into the right-hand sides of (p8|)-(^3|). This could be done 
without assuming that the magnitude of epjv is a first order term because the IPN effects have already been taken 
into account using the post-Keplerian parameterization instead of the Keplerian one. 
In celestial mechanics, the parameter T, is often replaced by (e.g. [|l]j6[0|): 

(j = -nT (34) 

which is just the constant related to the mean anomaly M ^ nt + a, i.e. the angle which the radius vector would 
have described if it had been moving uniformly with the average rate n. 

This newly defined orbital element allows one to see that all Lagrangian brackets are not changed except of this 
one: 



a jGm ( ^ Gm 
2c 



K«] = :772\/— (1 + —) (35) 



which takes the place of (P4). Using this new set of elements the IPN planetary equations become: 
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Usually the disturbing function Tii depends upon a only through the phase M = nt + a, therefore it is suitable 
to express Tii in terms of M . Furthermore TLi depends upon a both explicitly and implicitly through M (in fact 
n = n(a)), so that the derivative of the disturbing function taken with respect to a, can be written as follows: 

dHi f dHi\ , dHi dn 



da \ da J dM da ' 

the first term on the right hand side is the derivative of TLi with respect to a while keeping M (and then n) as a 
constant, and the second term takes into account the dependence of Tii upon a through n. From the above expression 
and using eq. (|3^ we get: 
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Now, since — eq. ( |3(^ ) becomes 

da _ 2c^a^ iGm Gm\ dHi 
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while the time derivative of M can be written 
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Eqs. ( [HD and (|4|) take the place of ( p6D and ( p^ ) when M is used as orbital parameter (instead of ct or T). 

V. APPLICATION 

As an example of application we consider the interaction of a close binary system with a third body. For the sake of 
simplicity we assume the third body at a distance Rb » a from the binary centre of mass with mass 7713 >> m. In 
this way the resulting motion could be safely approximated by the perturbed IPN motion of the binary system around 
its centre of mass which, in turns, has Keplerian orbit around TO3, considered motionless. The time evolution of 
rs is therefore known, this latter not being a variable of the problem. The discussion made at the beginning of section 
III showed that, in order to achieve a consistent picture, the effect of the perturbation must be greater than that of 
the 2PN one. The most interesting case is when the effect of the perturbation is of the same order of magnitude as 
the IPN effect; as it will be seen later this occurs when 



ma Gm^ / M^ 



where the second relation holds true if m ~ Mq and a ~ 2 i?© . 

Therefore, in the framework of previous assumptions, the perturbation to Hamiltonian (0) is given by 
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where Rb is the radius vector from my, to the IPN centre of mass of the binary system. We assume that the binary 
centre of mass performs a circular orbit around ms, which is coplanar with the orbital plane of the binary system. 
Hence we set i = 0, = and Rb = AB{cosijjyt, sincjat, 0), where 0^3 = ^J G ni^/A'^g. The Hamiltonian perturbation 
thus becomes 
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ni(a,e,M,uj,t) = - a^(l - ecosTy)^ [1 + 3 cos 2 (wg t - w ~ 77)] (48) 

4 A'^ 

where ry = 7/(a, e, M) and ry = 77(0, 77(0, e, M), e) (see eqs. (^) and (p^), together with the definition of M after eq. (|3^)) 
and 
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From eqs. ( p^ and (|l6|) one also has dn/da ~ —n/a, dn/de — 2eK/(l — e^), de^/da — {e — e^)/a, de^/de — es/e. 
In order to obtain the equations providing the temporal derivatives of the orbital elements, it suffices to calculate the 
partial derivatives of Tii; in this way the problem of motion is solved. In order to to better understand the behaviour 
of the solution let us suppose e = 0. In this case the Hamiltonian perturbation greatly simplifies 

Hi {a,M,Lu,t) = - ^„ [1 + 3 cos 2 (tj3 i - - (1 + k)M)] (50) 
4c^ Ag 

Finally the equations for the orbital elements a and M are: 
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while 

^ = 0; ^=0; ^ = 0; #=0 (53) 
dt dt dt dt ^ ^ 

that is e, uj, and i are constants of the motion. By substituting the unperturbed values of the orbital elements in the 
right-hand side of eqs. (|5l|)-(p2|) the first-order time-variation of the parameters are achieved. As far as relation ( |4^ ) 
is fulfilled, being nis a'^ /Jm ^ Gm/{c^ a), the solution can be written as 

a{t) = 00(1+ ^ ™^ "3° [cos 2 (w3 i - w - (1 + k)Mo) - cos 2uj]] (54) 

M{t) = (1- ^^^1 Mo + [sin 2 {uJ3t-uj-il + k)Mo) - sin2u;] (55) 

\ m Ag J 2 m Ag 

where ag and Mq are the unperturbed values. 

Accurate timing of the binary system could then provide the orbital profile through which (via eqs. (|5^)-(|55|)) 77^3 
and Ab are evaluated. This example has shown how the method proposed in this paper allowed one to achieve an 
already existing result in a relatively simple way. 



VI. CONCLUSIONS 



We have developed a Post-Newtonian extension of the planetary Lagrangian equations in order to describe the 
motion of a relativistic binary system under the influence of a generic external perturbation. An important class of 
phenomena that could be treated is the perturbation caused by external bodies to a binary system. We also have 
provided an example featuring the capability of the method. 

Our approach is a method to solve in a perturbative way the equations governing the variation of orbital parameters. 
The theory only assumes the smallness of the external perturbation while relativistic effects are already included in 
the unperturbed problem. This is the major advantage of our approach over classical Lagrangian methods. In the 
present paper the problem is solved at IPN level, but this approach should also be suitable to 2PN or higher-accuracy 
orders in the relativistic expansion. 
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